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We present the first comprehensive study of hybrid baryons using lattice QCD methods. Using a 
large basis of composite QCD interpolating fields we extract an extensive spectrum of baryon states 
and isolate those of hybrid character using their relatively large overlap onto operators which sample 
gluonic excitations. 

We consider the spectrum of Nucleon and Delta states at several quark masses finding a set of 
positive parity hybrid baryons with quantum numbers N^^+i ^1/2+1 ^3/2+1 ^3/2+1 -^5/2+1 anf l 
A 1//2 +, A 3/ / 2 + at an energy scale above the first band of 'conventional' excited positive parity 
baryons. This pattern of states is compatible with a color octet gluonic excitation having J p — 1 + 
as previously reported in the hybrid meson sector and with a comparable energy scale for the 
excitation, suggesting a common bound-state construction for hybrid mesons and baryons. 
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I. INTRODUCTION 

While QCD is the accepted underlying theory of 
hadronic physics, with hadrons being viewed as bound 
states of strongly coupled quarks and gluons, the role 
of excited gluonic fields in determining the spectrum 
of mesons and baryons remains unclear. QCD in the 
strongly coupled low-energy regime gives rise to a num- 
ber of interesting phenomena that are exhibited in the 
spectrum of mesons and baryons, such as the sponta- 
neous breaking of chiral symmetry yielding effective 'con- 
stituent' quark degrees-of- freedom. These quasi-particles 
have the quantum numbers of quarks, and there are clear 
signs of a spectrum of excitations corresponding to their 
relative motion. It seems odd then, considering that the 
gluonic field is strongly coupled to itself and to quarks, 
that there are no unambiguous signals for states featuring 
a gluonic excitation. 

The simplest place to look for gluonic excitations 
would seem to be in the isoscalar meson spectrum, where 
states built purely out of glue, the 'glueballs', could be 
present. This has proven to be very difficult in prac- 
tice [T] with suggestions that qq excitations mix strongly 
with glueball basis states. Another target is elsewhere in 
the spectrum of mesons, where states built from qq aug- 
mented by a gluonic excitation can have J quantum 
numbers not available to a pure qq state. These 'exotic 
hybrid mesons' remain our best hope of a 'smoking gun' 
experimental signature, but in practice the current ex- 
perimental situation is at best confused (see the review 
in [2]). Theoretical understanding of hybrid meson states 
in QCD recently took a significant step forward with the 
application of lattice methods to the problem. In [3H6] 
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the spectrum of mesons extracted from lattice QCD com- 
putations was interpreted in terms of 'conventional' qq 
states supplemented with a spectrum of hybrid mesons, 
some with exotic J c and some with conventional J PC . 
The degeneracy pattern of these states and the form of 
the composite QCD operators used to interpolate them 
from the vacuum strongly suggested that the gluonic ex- 
citation present in hybrid mesons is of chromomagnetic 
character. 

Hybrid baryons have not attracted the same attention 
as hybrid mesons principally because they lack manifest 
'exotic' character. All J p values can be populated by 
states constructed from three quarks having excitation 
in orbital angular momentum, so that hybrids can only 
appear in terms of over-population with respect to some 
model of qqq excitations. While a dynamical model may 
suggest peculiar decay characteristics for hybrid baryon 
states that might help to identify them, there is not the 
'smoking gun' quantum number signal present in the me- 
son sector. Furthermore the current experimental situa- 
tion, in which far fewer baryon resonances are observed 
than are expected in qqq bound-state models, does not 
encourage adding additional gluonic degrees-of-freedom 
to the bound-state system. 

Despite there being, as yet, no simple way to extract 
them experimentally, there is significant theoretical in- 
sight to be gained from studying hybrid baryons, par- 
ticularly in a framework that can simultaneously calcu- 
late the properties of hybrid mesons. We might imagine, 
given the different number of quarks, and the correspond- 
ing differing distributions of color sources, that the spec- 
trum of gluonic excitations in a baryon could be different 
to that in a meson. We will explore this possibility in 
this paper. 

The spectrum of hybrid baryons has previously been 
considered in a number of QCD-motivated models. In 
the bag model [7], in which quark and gluon fields are 
confined within a cavity with the fields satisfying ap- 
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propriate boundary conditions at the wall of the cavity, 
a TE-mode gluonic field (transforming as J PC = l" 1 ) 
combines with three quarks in an overall color octet to 
produce a lightest set of hybrid baryons : 

Nj./2+i N-L/2+, N 3 / 2 +, N 3 / 2 +, N 5 / 2 + 
A-1/2+, ^3/2+- 

The overall mass scale of these states is somewhat plas- 
tic, depending upon how one chooses to set parameters 
within the model - in [7|, the first excited pion resonance, 
the 7r(1300), was assumed to be a hybrid meson which 
then sets the scale for Nucleon hybrids starting as light 
as 1.6 GeV. A more plausible modern candidate for a 
hybrid pion would be the second excited pion resonance, 
the 7T (1800), which would place the hybrid Nucleons sig- 
nificantly heavier, somewhere above 2 GeV. 

The QCD sum-rule study in [8j considered the possibil- 
ity of a hybrid baryon as the first-excited state above the 
nucleon with Ni/ 2 + quantum numbers concluding that 
such a state would lie close to 1.5 GeV. 

An implementation of the flux-tube model EH], in 
which quarks sit on the ends of strings, providing a lin- 
ear confining potential, suggests a set of lightest hybrid 
baryons |11) : 

N\/2+, Nl/2+, A3/2+, N 3 / 2 + 
A1/2+, A3/2+, A 5 / 2 +- 

This differs from the bag model spectrum by having the 
J p = 5/2+ state with isospin-3/2 rather than isospin- 
1/2. The mass scale of these states is found to be around 
1.9 GeV. The flux-tube model is an example of a frame- 
work where we would not necessarily expect there to be 
a common energy scale for the gluonic excitation in a 
baryon and in a meson. A meson contains a single flux- 
tube between the quark and anti-quark which in a hybrid 
is excited in transverse oscillation. In a baryon however 
there are three tubes which either meet at a junction or 
form a triangle. There would seem to be no particular 
reason the excitations of this more complicated system 
need be the same as a single tube. 

As mentioned earlier, a major challenge for the ex- 
perimental isolation of hybrid baryons is that they have 
non-exotic quantum numbers and hence appear embed- 
ded in a spectrum of conventional baryon states. In fact, 
worse than this there is, a priori, nothing within QCD 
to prevent hybrids mixing arbitrarily strongly with con- 
ventional qqq states of the same J p to produced mass 
eigenstates which are neither one nor the other. This 
complication will also appear in any sufficiently sophisti- 
cated theoretical approach. Whether QCD in fact mani- 
fests strong mixing will need to be determined in explicit 
calculation. 

In this paper we will extract a spectrum of QCD eigen- 
states from baryon correlators computed using discreti- 
sation of QCD on a finite lattice. This technique is an 
ab-initio non-perturbative method starting from a gluon 



and quark lagrangian that allows us to compute the spec- 
trum directly from a controlled approximation to QCD. 
In order to explore the spectrum of states with a given set 
of quantum numbers it has proved effective [SHU |T2T - fri] 
to construct a large basis of hadron interpolating fields, 
featuring combinations of quark and gluon fields having 
the desired external quantum numbers. A matrix of two- 
point correlation functions can then be evaluated and by 
diagonalisation, a variational best estimate of the spec- 
trum of states determined. The eigenvectors obtained 
in this procedure indicate the optimum linear combina- 
tion of interpolating fields for each state in the spectrum. 
Quark-gluon bound-state structure interpretational in- 
formation follows if a particular state has a large over- 
lap with certain characteristic interpolating fields, and 
large overlap onto multiple interpolating fields of differ- 
ent characteristic structure might suggest large mixing of 
basis states within QCD. 

Our approach will mirror that used to determine a 
spectrum of hybrid mesons in [H 0J [6] , where we found 
that the lightest supermultiplct of hybrid mesons con- 
sists of an exotic J PC = 1 h state and three non-exotic 
states with J PC = '", 2 1 which are embedded 
in a spectrum of qq excitations. The character of gluonic 
excitation that appears to give rise to this spectrum is 
chromomagnetic, having J PC = l" 1 . Heavier positive 
parity hybrid mesons were also identified that appear to 
correspond to quark orbital excitation on top of the same 
chromomagnetic gluonic excitation with states featuring 
any other type of gluonic excitation apparently signifi- 
cantly heavier. 

The corresponding procedure for baryons will be fol- 
lowed in this paper: using a large basis of interpolating 
fields, including a number which rely upon the presence 
of a non-trivial gluonic field, we will extract a spectrum of 
baryon states. By considering the relative size of overlap 
of the extracted states onto interpolators of characteris- 
tically qqq and hybrid form we will identify those which 
we believe to be hybrid baryons. 

As with all lattice QCD calculations, we are somewhat 
limited in scope by computational restrictions which pre- 
vent us from calculating with quarks whose masses are 
as light as those in the real world. In place of this we 
will compute with a range of heavier quark masses and 
study the trend as the quark mass is reduced. Since the 
computational cost increases with the volume of space- 
time considered we limit ourselves to a small region a few 
fermi across discretised with grid points a fraction of a 
fermi apart. 

In this work we will report on the first systematic stud- 
ies of hybrid baryons using lattice QCD methods. We 
will work with three flavours of dynamical quarks - the 
strange quark being fixed at its physical mass and the 
degenerate up and down quarks being varied down from 
the strange quark mass (where we would have an exact 
SU(3)f symmetry) to a value where the pion mass is 
400 McV. This is an extension of the work presented in 
[15] which computed the baryon spectrum on the same 
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lattices, but excluded baryon interpolators featuring a 
gluonic excitation and which made no observations re- 
garding hybrid baryons. 



II. BARYON INTERPOLATORS & THEIR 
INTERPRETATION 

The construction of a basis of baryon operators re- 
specting the required permutational symmetry of three 
quark fields and featuring up to two gauge-covariant 
derivatives was described in [15]. This leads to a rather 
large basis of operators that covers all spins and both 
parities up to J p = | . For lower spins there are suffi- 
ciently many operators that variational analysis can lead 
to extraction of a significant number of excited states. 

The baryon interpolators used can be expressed in a 
compact notation which exposes the permutational sym- 
metries, 5If,s,d = S, M, A, 

(B Tf ®{S p *)l®Df TD )\ (1) 

where the three factors describe the flavor, Dirac spin and 
derivative structure of the interpolator. In this report we 
will consider only baryons of flavor B-£ f = Nm, As, the 
strangeness-zero members of 8^ , 10 f representations, 
over a range of quark masses. The only exception to this 
is the special case of three degenerate flavors of quark, 
i.e. having exact SU (3) flavor symmetry, where in ad- 



dition to the octet and decuplet states, we'll extract a 
spectrum of B-£ F = Aa states, the only member of the 
1^ representation. 

The Dirac spin factor, (S^) , represents the possi- 
ble combinations of upper and lower components of three 
Dirac spinors transforming with angular momentum S 
and parity P$ and having permutational symmetry, Us- 
For some {S Ps )y. s there are multiple possible construc- 
tions and they are labelled by the superscript n. 

The derivative factor, D^t d > indicates the number of 
gauge-covariant derivatives applied to the quark fields 
(d = 0, 1, 2) and that they are combined to give objects 
transforming under rotations with angular momentum L 
and with permutational symmetry Zq. 

With a single derivative (d = 1), excluding a symmet- 
ric total derivative, only mixed symmetrjjj possibilities 
transforming as L = 1 arise: 

f D [1] - 1 (W {3) - n (1) - n (2) ) 

n [l] I U MS,m - V> u m u rn U m ) 

U L=\M ■ \ n [l] _ W n (l) _ n (2)\ 
L MA.in — ^2 V m m J 

where e.g. fl™ = e(m) ■ D^> is the gauge-covariant 
derivative, in a circular polarisation basis, acting on the 
first quark field. 

At the two derivative level (d = 2) definite permuta- 
tion symmetry is implemented by the appropriate lin- 
ear combination of products of -D^g ma , while the rota- 
tional transformational properties are set using an SO (3) 
Clebsch-Gordan coupling forl<g)l->L = 0ffil®2: 



D 



[2] 
L,M 



MS,m' 



MA,m MA,m' / ' 



4 2 ; s = (lm;lm'\LM)^=(D l ^ m D 

4 2 !ms = (1m; lm'\LM)^(-D^ Stm D^ s ^ + D^ A m D^ A m ,), 
D [ l] M a = (lm;lm f \LM)^(D^ Stm D^ 

D^ A = (lm;lm'\LM)^(D^ 5<m D 



i] _ D [l] 

MA,?-' 1J ' 



D 



[i] 



MA,m MS,m 



(2) 
(3) 
(4) 



In practice, the interpolators ([I]) correspond to sums 
of terms of generic structure 

£abc (>4(l± 7 >) a (> 2 i(l±7>) 6 (-D" 3 |(l±7 )^ 

where each quark field is projected into either upper or 
lower Dirac components and is acted upon by either zero, 
one or two gauge-covariant derivatives. The spin and fla- 
vor couplings are suppressed here (details of the construc- 
tion can be found in |15j). The required antisymmetry 



MS, MA indicates symmetric or antisymmetric behavior in 1 o 2 



in color, producing overall color-singlet operators, is im- 
plemented using the Levi-Civita symbol. Given this it is 
clear that the flavor, spin and spatial symmetry repre- 
sentations, Zf,5Is,5Id must be combined to give overall 
symmetric (S) combinations in order for the interpolator 
to obey Fermi antisymmetry. 



A. "Conventional" qqq interpretation of operators 

A notable subset of the interpolators constructed are 
those featuring only upper components of Dirac spinors, 
which we refer to as "non-relativistic" . As shown in Table 
IV of [TS] , these operators can be simply classified accord- 
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ing to an SU (3) f®SU{2)s®0{3>) symmetry analogous to 
what one might have in a (non-relativistic) three-quark 
bound-state model. In this assignment, the derivatives 
play the role of orbital angular momentum between the 
quarks (the 0(3) symmetry in most bound-state models). 

As a concrete example consider the symmetric con- 
struction of two derivatives defined in equation (J2J) , D 
which can be explicitly written 



2] 

L,S> 



D 



[2] 



3^2 



(lm; lm'\LM) 



l( £>« D (1) + D (2) + D (3) ] 

m m' 1 m m' J 



(5) 



(3) (2) 

where expressed more fully a term D m D m , would have 
the color-spatial structure e a b c tp a (D m /ip) b (D m ip) c . A 

[21 

simple model interpretation of the action of D L S fol- 
lows if we treat the gauge-covariant derivatives as or- 
dinary, commuting, derivative operators (neglecting the 
presence of the gauge-field) and use Jacobi co-ordinates 
(p = -^(n - r 2 ), A = -L(r*! + r 2 - 2f 3 )) for the quark 
positions. In that case we find that 

d l]s ~ ( lm ! lrn'\LM) (p p ) m (Pp)m' + (p\)m(p\) m > 



[p 2 p Yi I (p P )YS{px)+plYt I {px)YS(p p ) 



^L,cvcn ; 



so that if L — 2, we could interpret this as the symmetric 
combination of a D-wave on the p- "oscillator" coupled to 
an 5- wave on the A- "oscillator" and a D-wave on the A- 
"oscillator" coupled to an S- wave on the p- "oscillator" . 
This kind of operator, coupled to "non-relativistic" Dirac 
spinors can interpolate from the vacuum exactly the kind 
of structures used in qqq "quark models" 16-19J. For 

3 + 



[2] 

L=2,S 



example, the interpolator (jVm ® (^ + ) M ® D 

would be identified with the model state N 2 D$ | + . Table 
IV in 15J presents qqq model interpretations for "non- 
relativistic" interpolators. 



B. "Hybrid" qqqG interpretation of operators 

The subset of operators which we identify as essentially 
hybrid in nature appears at two derivatives and follows 
from combinations which correspond to commutation of 
two gauge-covariant derivatives acting on the same quark 
field. Such operators are zero in a theory without gauge- 
fields and correspond to the chromomagnetic components 
of the gluonic field-strength tensor. This is the central 



observation that leads us to identify these operators as 
being 'hybrid' - unless an eigenstate contains a gluonic 
field in a non-trivial configuration it should have no over- 
lap onto an operator of this form. 

In the derivative constructions described above 
such entries appear via terms of the form 

(lmjlm'llAOD^D^ <x e(M) i e iife [£» i)J D fc ] cx B M , 
where the chromomagnetic field in a circular polarisa- 
tion basis, Bm = e{M) ■ B is related to the commutator 
of two gauge-covariant derivatives, Bk = o e ijk[Di, Dj]. 
As an example of how such an object enters into our 
operator basis, consider the projection of equations ^ 
into L = 1: 



D 



[2] 

L = 1,MS 



3V2 



(lm; lm'| 1M) 



2 ( n(i) D w + D m D m _ 2D ( 3 ) £,0) 

' D {3 } + D {1 }) 

1 - m 1 - vr\ ' 1 777 TT-i 1 



-2(dMd® + d®d<3 



D 



L = 1,MA 



V6 



(lm; lm'jlM) 

( -D^D (1) ,+D^D (2 ] 
+ (D^D {3) , +D^D {1) , 
- (D^D {3 1 +D ( ^D {2) , 



In each case there are terms which are symmetric in 
m m' while (lm;lm'|lAf^ is antisymmetric - since 
derivative operators acting on different quarks commute, 
these terms are zero. The remaining terms, in which both 
derivatives act on the same quark field, give color-spatial 
structures proportional to 

MS ~ e abc ({B M i>)aipb^c + ipa(B M ip)bipc - 2ip a ip b (B M ip)c) , 
MA ~ e abc (-(B M ip)a^b^c + ipa{B M 4>)bipc) ■ 

The form of the three-quark color coupling can be ex- 
posed by expressing the chromomagnetic field in an ad- 
joint basis, Bm = B M t A , so that 



MS 



^abc {-~{t A ^)a^b1pc + ?Pa{t A l(>)b4>c) , 



which indicates that the three quarks are coupled to a 
color-octet as either (3 ® 3 -> 3) ® 3 ->• 8 — MS or 
(3®3->6)®3->8~ MA . The chromomagetic field (a 
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L = 0,2 


L = 1 




(<7W)ic 





M 




[(9W)8 c G 8c ] lc 







(W9)ic 



TABLE I. Summary of simplest bound-state interpretations 
of two-derivative structures. In practice the derivatives are 
actually implemented as finite differences through gauge- 
covariant parallel transport (links) and this gives rise to cor- 
rections to the above assignments that are suppressed by pow- 
ers of the lattice spacing. 



color octet) is coupled in to give an overall color singlet 
baryon operator, (il'4'i') A B^f. These are the operators 
which we identify with hybrid baryons. 

At first glance it would appear that Df^ =1 s , as defined 
in equation ^ and written out in full in equation ([5| 
should also have hybrid character owing to the presence 
of terms featuring (lm; lm'\lM)Dm^ D^) . In fact this 

[21 

operator is zero - D L=1 s gives rise to a three-quark color 
structure 

which cannot be a color-octet as the antisymmetric 
combination of three 3 C representations gives a color- 
singlet. We do not include any operators constructed 
using D^ =1 s in our basis. 




where ^g, ipe are the MA, MS color combinations (cor- 
responding to the MS, MA derivative combinations) and 
0, x are the flavor, spin structures presented in Appendix 
A of [15] ■ This set of operators corresponds to a [70, 1 + ] 
in the SU(6) fs assignment scheme, a representation not 
present for conventional three-quark baryons. 

[2] 

Coupling the hybrid derivative structures (D L=1 M ) to 



We mention in passing that the antisymmetric com- 
bination of two derivatives defined in equation Q, has 
structure 

D f = ^ ([D^\D^] + [DW,DW] + [dW,dW]) , 

which does not feature the commutator of two gauge- 
covariant derivatives acting on the same quark field. The 
antisymmetric nature is such that only with L = 1 do we 
get non-zero operators but these are not of hybrid na- 
ture; in the Jacobi basis presented in the previous section 
they transform like (lm; lm' \1M} (p p ) m (p\) m i , i.e. they 
can be viewed as a P-wave on each "oscillator" in the 
language of the constituent quark model, giving states 
labeled 2 ' 4 P A - 

The interpretations of the constructions are sum- 
marised in Table [TJ These operators, excluding those iden- 
tified above as being of hybrid character, were used in [T5J 
to extract a spectrum of baryons. In the next section 
of this paper we will present spectrum results including 
the hybrid operators, identifying a number of states with 
large overlap onto these gluonic operators that we inter- 
pret as being hybrid baryons. 

[21 

We have determined that operators featuring D L=1 M 
can be associated with hybrid structure - now we can 
ask what flavor-spin constructions can be combined with 

[21 

D L=1 M to give operators which respect Fermi antisym- 
metry? As is the case with qqq interpretations, the easiest 
hybrid interpolators to interpret phenomenologically are 
those with "non-relativistic" Dirac spin structure. The 
explicit flavor-spin-color structure of these operators, fea- 

[21 

turing D L=X M , is as follows: 



~ I [^a (0maXma - </»msXms) - "06 (</>maXms + 0msXma)] 
= 73 [4>3 Ox)ms - Ox)ma] = W4 

~ 75XS (V>30MS - -06</>Ma) = XS (lp<t>) A = [^0X] A 

~ TS^s (1P3XMS ~ -06Xma) = 0s ("0x)a = WxIa ( 6 ) 

I 

Dirac spin lower components gives a somewhat larger set 
of operators, including some with negative parity, as indi- 
cated in Table [nj although their model interpretation is 
somewhat more involved and will not be described here. 

In order to account for the reduced (cubic) rotational 
symmetry of the lattice, we subduce all interpolators into 
irreducible representations of the cubic group and form 
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positive parity 


negative parity 


N(} + 


) 32{4}[8](2) 


H¥ 


h ) 16{2}[4](1) 




32{2}[8](0) 




16{1}[4](0) 


N(l + 


) 36{5}[8](2) 


A(|- 


h ) 19{3}[4](1) 


AT(f-) 


36{2}[8](0) 


A(F) 


19{1}[4](0) 


N(l + 


) 19{3}[3](1) 


A(f- 


h ) 9{2}[1](0) 


iv(f-) 


19{1}[3](0) 


a(F) 


9{0}[1](0) 


N(l + 


) 4{1}[0](0) 




h ) 3{1}[0](0) 




4{0}[0](0) 


A(F) 


3{0}[0](0) 



TABLE II. Number of operators constructed in each quantum number channel: total{non-relativistic conventional} [hybrid] (non- 
relativistic hybrid) 



correlators according to that symmetry. In [15] we show 
that the spectrum of states extracted from variational 
analysis of these correlators can be spin-identified with 
confidence and herein we will present the spectra labelled 
according to the continuum spin. 

III. BARYON SPECTRUM 

Mass spectra and matrix-elements follow from vari- 
ational analysis of a matrix of correlators, Cij{t) = 
(0|Oi(t)Oj(0)|0>, which can be shown to correspond to 
solution of the generalised eigenvalue problem 

C(t)v a = X n (t,t Q )C(t )v n . 

Each state |n) that can be interpolated by some lin- 
ear combination of 0\ has associated with it a "prin- 
cipal correlator", X n (t,t ) that at large times behaves 
like e~ E "^ to \ yielding the energy of the state, and an 
eigenvector v n which indicates the optimum interpolator 
for the state, J2i V V®\- Details of our procedure for effi- 
ciently and accurately implementing such a solution can 
be found in [H EQ]. 

To the extent that one has a bound-state interpretation 
of the interpolator Oi, built from quark and gluon fields, 
the relative sizes of the matrix elements (n|cj(0)|0) can 
be used to interpret the bound-state structure of state 
|n). These matrix elements are trivially related to the 
eigenvectors of the generalised eigenvalue problem solved 
above. Such a procedure was followed in [J], leading to 
identification of the low-lying meson spectrum as being 
rather like the constituent quark model qq spectrum, sup- 
plemented with a spectrum of hybrid mesons in which 
a chromomagnetic gluonic field was dominant. We will 
take the same approach here in the baryon spectrum in 
order to isolate the role of hybrid baryon basis states in 
the spectrum. 

We performed calculations on lattices of size 16 3 x 128, 
corresponding to a physical spatial volume of approxi- 
mately (2.0 fm) 3 , with three flavors of dynamical quarks. 
One quark has mass corresponding to the strange quark 
and the other two correspond to degenerate u and d 
quarks. One calculation had all three quarks degenerate 
at the strange quark mass, corresponding to 702 MeV 
octet pseudoscalar bosons, and the remaining two cal- 
culations had lighter u, d quarks giving pion masses of 
524 and 396 MeV. Details of these dynamical anisotropic 



clover configurations can be found in [STJ [23] . The two- 
point correlators used in the variational analysis were 
computed via distillation [23 in which all quark fields are 
smeared over space by an operator □ = Yl n con ~ 
structed using the lowest N = 56 eigenvectors of the 
three-dimensional gauge-covariant laplacian (— D 2 £ n — 
Ari^n)- All gauge-links entering in the operator construc- 
tions are stout-smeared [2"4"j . 

To facilitate comparisons of the spectrum at different 
quark masses, the ratio of hadron masses with the f2 
baryon mass obtained on the same ensemble is used to 
remove the explicit scale dependence, following Ref. |22) . 

A. Nucleons & Deltas at m v = 524 MeV 

Figure [l] shows the spectrum of positive parity Nucleon 
and Delta states extracted on the 524 MeV lattice along 
with histograms showing the relative sizes of matrix ele- 
ments (n|Oj(0)|0) for a set of "non-relativistic" interpo- 
lators characteristic of qqq and hybrid basis states. The 
histograms are normalised such that the largest overlap 
onto a given operator across the entire extracted spec- 
trum has unit value (the tallest bars in the figure are of 
unit value). 

Focussing first on the J p = | nucleon states we ob- 
serve a light ground state nucleon having dominamj^Jover- 
lap onto a "non-relativistic" interpolator (]Vm (g> (| ) M ® 

D L _ s ) 2 with qqq structure 2 5*s. The band of four ex- 

cited ^ states between 2.0 and 2.5 GeV can be char- 
acterised as being dominantly admixtures of the qqq ba- 
sis 2 Ss, 2 Sm, 4 -Dm, 2 Pa, and while this basis is clearly not 
manifested diagonally in the spectrum, the mixing is 
modest. Notably, the extracted spectrum lacks a light 
first-excited state that we might associate with the broad 
enhancement observed experimentally and known as the 
"Roper" resonance. We will return to this later when we 
discuss the role of hadronic decays on the spectrum. 

Slightly heavier than the first-excited band of four 

1 + 

2 states we observe two states with dominant over- 
lap onto non-relativistic interpolators we have identified 



2 The observed sub-dominant overlap of the ground state nucleon 
onto operators featuring an excited gluonic field is also required 
to describe QCD sum-rules of the N^i 2 + channel [8] 
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FIG. 1. Extracted spectrum of Nucleon and Deltas states by J p at a pion mass of 524 MeV. Rectangle height indicates the 
statistical uncertainty in mass determination. Histograms indicate the relative size of matrix elements (n|Oj(0)|0) for a "non- 
relativistic" subset of the operators used - normalisation is as described in the text with the lighter area at the head of each bar 
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and 4 hyb = {Nu ® (| + )s ® ^l=i m)- Asterisks indicate states having dominant overlap onto hybrid operators. 



as being of hybrid nature, (N M <8> (5 ) M ® D\L X M ) 2 

1 + 

and (A^m f? 1 (| + )s ® -^lLi m) 2 • They appear to be 
largely admixtures of two basis states with three-quark 
spin S = 1/2, 3/2 that we label 2 hyb, 4 hyb. We mention 



in passing that in [15], the spectrum extracted without 
hybrid interpolators showed one badly-determined state 
in this same mass region, indicating the importance of 
including interpolators with good overlap onto all types 
of state present in the spectrum. 
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The next nucleon state resolved with J p = A is some- 
what heavier than the hybrids and has no obvious inter- 
pretation in terms of overlaps with the "non-relativistic" 
operators used (with up to two gauge- covariant deriva- 
tives) . Still heavier states are present with these quantum 
numberifO 

In the J p = | nucleon channel, we see a lightest 
band of states roughly degenerate with the first-excited 
band in ^ + . This band features five states observed to 
be dominated by the qqq basis 4 5m, 2 Ds, 4 L>m, 2 Dm, 2 Pa, 
with the state dominated by 2 P& being heaviest as it was 
for J p = \ + ■ As in the nucleon ^ + case we observe 
a pair of states somewhat above this band which have 
dominant overlap onto 2 hyb, hyb. A few hundred MeV 
heavier there are more states of unidentified structure. 

In the J p = | nucleon channel there are three states 
in the "first-excited" mass region, apparently constructed 
from superpositions of 4 D$, 2 L>m, 4 L>m- Somewhat heav- 
ier is a single state having large overlap onto a hybrid 
operator of character 4 hyb. Unidentified | + states are 
present a little heavier than the hybrid. 

In the J p = ^ + nucleon channel, there is a single 
state in the "first-excited" region having dominant over- 
lap onto 4 Dm- Heavier states of unidentified structure 
are present above 3 GeV, with no identified hybrid states. 
We remind the reader that there are no "non-relativistic" 
hybrid interpolators having J p = | constructed at 
level. 

In the Delta sector we again see a band of "first- 
excited" states whose population and overlaps are in ac- 
cord with expectations of a qqq model. States with dom- 
inant overlap onto our hybrid interpolators appear above 
that band, but below other unidentified states. The hy- 
brids appear in the J p = ^ + , | + channels and are of 
2 hyb character. 

The spectrum of negative parity nucleons and Deltas 
below 2 GeV is almost identical to that presented in [TS] , 
obtained without inclusion of any hybrid interpolators, 
and agrees with qqq model state counting. Negative par- 
ity states having significant overlap onto the hybrid op- 
erators in the basis first appear above 3 GeV, somewhat 
heavier than the positive parity hybrid states. 

B. If, 8f, 10f at an SU(3)f symmetric point 

The flavor singlet construction in equation (J6J), Ai,a <8> 

(I + )m ® ^l=i m> i s mos t simply explored at an SU(3)f 
point where the states interpolated cannot mix with 
Lambda states in an 8^ representation. In Figure [2] we 
show the spectrum extracted with all three quark masses 



set to the strange quark mass - here the lightest pseu- 
doscalar meson has a mass of 702 MeV. We show the 
mass spectrum for 8^, 10^, lp representations. 

We observe that the 8^, 10j? spectra strongly resem- 
ble the Nucleon, Delta spectra extracted at the 524 MeV 
pion mass shown in figure [T] The 1^ spectrum consists of 
the states expected in a qqq model, lying roughly degen- 
erate with the first excited band in the octet, decuplet 
spectrum, and in addition, two states of hybrid character 
having J p = ^ + , | + . As in the octet and decuplet, the 
hybrid states are slightly heavier than the qqq states. 



C. Quark mass dependence 

In order to determine if the hybrid spectrum might 
be strongly quark mass dependent we also computed the 
Nucleon and Delta spectrum with u, d quark masses such 
that the pion weighs 396 MeV. The results are presented 
in Figure [3j where we see that the gross structure of the 
spectrum is as it was at ra, = 524 MeV. 



IV. CONCLUSIONS 

From the extracted baryon spectra presented in the 
previous section we can draw a number of conclusions 
regarding the nature of hybrid baryons and more gener- 
ally gluonic excitations within QCD. The states identified 
have large overlap onto interpolating fields containing the 

[21 

operator D L=1 M which transforms like a chromomag- 
netic field (color octet, J PC = U Q The form of the 
operators suggest that within a hybrid baryon the three 
quarks are arranged in a color octet with the chromomag- 
netic gluonic excitation making the state an overall color 
singlet. The low- lying hybrid states overlap strongly onto 
the "non-relativistic" subset of our hybrid interpolators, 
those constructed using only upper components of Dirac 
spinors. We can interpret this as suggesting that the 
quarks within the lightest hybrid baryons are dominantly 
in S- waves. 

The particular set of flavor- J p states observed, 
Ni/ 2 +, N 1/2 +, N 3/2 +, N 3/2 +, N 5/2 +, plus A 1/2 +, A 3/2+ 
and A lil /2+ , A 13 / 2 + , suggests that all color singlet states 
are present that can be constructed from antisymmetric 
(qqq)8 c with quarks in an S'-wave coupled to a chromo- 
magnetic (8 C , 1+) gluonic excitation (see equations Q). 
Furthermore the fact that we are easily able to pick the 
hybrid states out of a dense spectrum of qqq states indi- 
cates that they tend to retain their character despite the 
possibility of mixing strongly with qqq states of the same 
J p . 



3 There are a total of 32 i states extracted, equal to the number 
of interpolators used, but the statistical precision and stability 4 see [6] for a simple model in which this operator interpolates a 
of extraction of the masses decreases as one goes up in energy quasi-gluon in a P-wave 
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FIG. 3. Extracted spectrum of Nucleon and Delta states at a 
pion mass of 396 MeV. Grey boxes are conventional qqq states 
and blue boxes are the states identified as hybrid baryons. 



The position of hybrid baryons within the spectrum of 
expected qqq states is now determined: the lightest states 
are of positive parity and lie slightly heavier than the 
first-excited positive parity qqq states. Negative parity 
hybrid baryons appear to be heavier than this. 

In [5] we presented the spectrum of hybrid mesons ex- 
tracted on the same lattices. The distribution of states 
across J PC (the supermultiplet structure) is also com- 
patible with a color octet chromomagnetic field coupled 
to quarks in a color octet (in this case (qq)g a ). This 
leads us to think that the gluonic excitation form may 
be common to hybrid mesons and baryons. To deter- 
mine whether the energy scale of the gluonic excitation is 
common we choose to plot the spectrum of hybrid mesons 
alongside the spectrum of hybrid baryons. Of course for 



this to be meaningful we must take some account of the 
differing number of quarks - in a constituent quark model 
we would subtract twice the constituent quark mass in 
mesons and three times the mass in baryons. With the 
lattice data we opt to subtract the p mass from the meson 
spectrum and the nucleon mass from the baryon spec- 
trum. This is presented in Figure [4] where it would ap- 
pear that there is a common energy scale, in the region 
of 1.3 GeV, for the lowest gluonic excitation in mesons 
and baryons. 

It is worth mentioning here that in the meson sector we 
used an operator basis featuring up to three derivatives. 
The low-lying spectrum of negative-parity hybrid mesons 
was largely insensitive to the presence or absence of the 
three-derivative operators. Conversely, accurate determi- 
nation of the higher-lying positive-parity hybrid mesons 
required inclusion of these operators. We expect that a 
similar situation would hold true in the baryon spectrum 
- where we to include a set of three-derivative opera- 
tors, the low-lying positive parity hybrid baryon spec- 
trum would be largely unchanged, but we would likely 
resolve more clearly a spectrum of higher-lying negative 
parity hybrid baryons. 

There are features of our extracted spectrum that we 
cannot explain without more detailed modelling, in par- 
ticular the nature of the dynamics giving rise to the fine 
structure in the spectrum is not known to us. It is inter- 
esting to note that the extracted structure is rather close 
to that presented in [7J following from perturbative cor- 
rections (quark-spin, gluon Compton etc ...) to the bag 
spectrum. One possible exercise using the spectrum ob- 
tained in this paper would be to explore the minimal set 
of model components required to describe the observed 
fine structure. 
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With optimum linear combinations of baryon interpo- 
lators determined in the variational solution of the matrix 
of two-point correlators, we can in the near future con- 
sider computations of radiative transition form-factors. 
Such calculations, which involve three-point hadron cor- 
relators, have been done for charmonium |25| . and for 
spin- 1/2 baryons, albeit with very small operator bases 
[2"fJl 127) . Radiative transition form- factors are of phe- 
nomenological interest because they can be measured 
in baryon electroproduction experiments [28^ where the 
photon Q 2 dependence and the relative size of multipole 
amplitudes might provide some insight into the internal 
structure of the excited state [29] . 

The physics of excited states as resonances, able to de- 
cay into meson-baryon final states is incomplete in this 
study. Along with larger than physical quark masses, this 
may be the cause of the lack of a light Roper candidate 
state within the calculation. In order to observe such 
decay physics in the extracted finite volume spectrum, 
which would correspond to additional volume-dependent 
energy levels, we must supplement our basis of inter- 
polating fields with some resembling multi-hadron con- 
structions (e.g. nN, ttttN). Such work is underway and 
may eventually lead to phenomenological decay filters for 
hybrid character that would be useful for experimental 
searches. 

A more complete understanding of hybrid baryons in 
QCD will of course require calculations at lighter quark 
masses, approaching the true physical values, coupled 
with larger lattice volumes to accommodate the increas- 
ing physical size of the bound states and such work is 
now warranted given the success of this first calculation. 



In summary we have extracted a spectrum of baryons 
from QCD computations that features 'expected' qqq 
states built from superpositions of flavor-spin representa- 
tions [56,0+], [70,0+], [70,2+], [56,2+], [20,1+] in posi- 
tive parity and [70, 1 _ ] in negative parity, supplemented 
with a set of positive parity hybrid baryons that form 
a [70, 1+] and which lie above the band of first-excited 
positive parity baryons. In light of this it seems unlikely 
that the Roper is dominantly of hybrid character as has 
been speculated in the past. The structure of the opera- 
tors that interpolate the hybrid states from the vacuum, 
along with the observation that the energy scale of glu- 
onic excitations appears to be common for mesons and 
baryons, provides evidence that the gluonic excitation 
sector in QCD may turn out to be relatively simple. We 
suggest that a chromomagnetic excitation (J PC = 1 H ) 
is lightest with an energy scale in the region of 1.3 GeV. 
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